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ABSTRACT 

We present the results of an archival search for Trans-neptunian objects (TNOs) in an ecliptic 
field observed with Subaru in 2002. The depth of the search allowed us to find 20 new TNOs with 
magnitudes between R — 24 and 27. We fit a double power law model to the data^ the most likely 



values for the bright and faint power law exponents are Q!i=0.73lQ Qg and = 



_Q 12 ^-nd the break magnitude is Req 



0.20t!^:i4; the differential 
This is the 



=25.0tg:«. 



number density at i? = 23 is (723=1. 
most precise measurement of the break in the TNO luminosity function to date. The break in the size 
distribution corresponds to a diameter of I? = 90 ± 30 km assuming a 4% albedo 
Subject headings: Kuiper Belt - Solar System: formation 



1. INTRODUCTION 

The TNO size distribution is the result of the forma- 
tion and coUisional history of its members. Models of 
the physical evolution assume the processes at play re- 
sult in a cha racteristic power la w size distribution. In 
the model of iPan &: Saril (j2005[ ) of strengthless TNOs, 
the collisional evolution of the objects leads to an evolv- 
ing size distribution that changes slope at the size at 
which the collisional lifetime of an object equals the age 
of the system. This analytical res ult confirms those ob- 
tained from numerical simulations (iKenvon fc Luiilll999t 
iDavis et al.lll999tlKenvon fc Bromlevll2004f) . The precise 
location of the break in the size distribution is an impor- 
tant measurement to link TNO formation and evolution 
models to the current TNO population fsee lKenvon et all 
[20081 for a review). 

The TNO luminosity function is defined by the dis- 
tance, albedo, and size distributions of TNOs, and it 
is readily observed. The distance distributio n exhibits 
a sharp edge at 50 AU ()Truiillo et"all l2001h that has 
been observed in other surveys (jCladman et all l2001t 
iFuentes fc Holmanll2008D . Though it is customary to use 
a 4 % albedo, correlations w ith size have been reported 
(see IStansberry et al.) 12008^ for a review) . After mak- 
ing sensible assumptions for the albedo and distances of 
the observed TNOs, the size distribution can be derived 
from the luminosit y function. S everal sea rches have been 
comp leted (see Be rnstein et al.l [2004: Fue ntes fc Holmanl 
20081 for a review), the deepest being the search by 



Bernstein et al.l (|2004l ) with the Hubble Space Telescope 



(HST). While previous results sampled the size distribu- 
tion up to magnitude R = 26 and were consistent with 
a pow er law luminosity function, the iBernstein et al.l 
(jlpoi) survey found a significant underabundance of ob- 
jects at i? ~ 28, compared to that predicted from extrap- 
olating the luminosity function at brighter magnitudes. 
They concluded a break in the luminosity function must 
occur between i? ~ 25 and 28. 

^ Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA; cfuentes@cfa.harvard.edu 

^ Current Address: Department of Astronomy, University of Cal- 
ifornia at Berkeley, Berkeley, CA 94720, USA 



Since IBernstein et al.l (|2004f ) announced their results, 
a number of ground-based surveys for fainter TNOs 
have been conducted in an effort to bridge the di- 
vide between the brighter (i?~23) TNOs found through 
wide -field searches su ch as the Deep Ecliptic Sur- 
vey (jElliot et al.l 120051) and the C anada-France Eclip- 
tic Plane Survey (iJones et al.ll2006l) and tho s e foun d by 
IBernstein et~aIT (|2004[ ). iFuentes fc Holmanl (|2008l ) dis- 
covered 82 TNOs in an archival search of Subaru data. 
This survey reached a limiting magnitude of i? = 25.7 
and successfully detected the break in the luminosity 
function. In the present work we build upon that ear- 
lier work, increasing the limiting magnitude of ground- 
based surveys to i? ~ 27 and further narrowing the 
gap between ground-based and space-based TNO sur- 
veys. Constraints on even fainter ob jects have been 
placed by stellar occultation surveys (|Bickerton et al.l 
I2OO8I: iZhang et al.i r20081 . 

We present the details of the data in the next section. 
§ [3] outlines the method used to analyze and search for 
TNOs in the data. The characterization of the search us- 
ing a synthetic population is presented in § [H In § [5] and 
§ [S] we present our results and discuss what they imply 
for the formation and evolution of the TNO population. 

2. DATA 

The apparent motion of solar system bodies limits the 
useful exposure time to the period required for an ob- 
ject to move a PSF width; longer exposures spread the 
signal along a trail and, correspondingly, increase the 
contribution of the sky background. To overcome this, 
we consider a series of short exposures and shift the suc- 
cessive images to compensate for the motion of the ob- 
ject. In this way the signal from the source can domi- 
nate the noise from the background. This "pencil beam" 
or "digital tracking" approach has been successfully ap- 
plied in searches for TNOs and outer planet satellites 
( Alle n et al.l 120011: iG ladman ct al. 2001; Hol man et all 
120041: iKavelaars et al..,2004: Eraser et al.,.2008i) . 

Usi ng Subaru's electronic archive SMOKA (jlchikawal 
l2002f) . we identified a data set well suited to such a 
search, a series of 148 consecutive 120 s exposures of 
a single, ecliptic field observed with Subaru on UT 2002 
September 2 over the course of 8 hours. These dat a were 
collected with Suprime-Cam (jMivazaki et al.|[2"002D . a 10- 
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CCD mosaic camera with a 34' x 27' field of view and an 
image scale of 0'.'202 pix~^. It is our understanding that 
these data were originally colle cted for the Subaru Mai n 
Belt Asteroid Survey (see Yos hida fc Nakamural l2007l ) . 
The field, R.A. 22:41:38, Dec -07:37:35, is less than 1° 
from the ecliptic and was observed near opposition using 
the "Cousins R" filter. Although brief intervals of clouds 
are apparent from the images, the night was mostly pho- 
tometric, with stable seeing of ~ 0'.'7. 

We trimmed, bias subtracted and flat divided the im- 
ages with calibrations obtained on the same night and 
the next one using standard IRAF'* routines. 

3. MOVING OBJECT DETECTION 

We determined an astrometric solution for each in- 
divid ual CCD using the 2MASS point source cata- 
log (Cutri et al.ll2003D . The RMS of the solution was 
~ 0'.'2, comparable to that of the catalog and much 
smaller than the typical seeing. We then registered every 
image to the same astrometric reference. 

Next, we inserted a population of synthetic objects 
that will be later used to determine the efhciency of our 
search (the details of this process are given in § |4]) . 

To detect TNOs moving at any physically plausible 
velocity, we defined a grid of rates, parallel and perpen- 
dicular to the ecliptic. A TNO at a typical distance of 
42 AU exhibits a parallactic motion of 3 " h~^, mov- 
ing about 120 pixels over the length of the observations. 
We searched rates from 0'.'7 to 5'.'1 h~^ for the paral- 
lel rate to cover the range between ~ 20 and 200 AU, 
and perpendicular rates in the range — 1'.'4 to 1'.'4 h~^. 
Furthermore, we restricted our attention to directions of 
apparent motion within 15° of the ecliptic, which nev- 
ertheless permits the detection of highly inclined TNOs. 
We searched these rates with a resolution of O'.'l h~^ 
along both axes (this resolution ensures that the signal 
in the first image can be aligned with that in the last for 
any object). 

Pr ior to combin ing the images, we used the ISIS pack- 
age (jAlardl l2000f ) to PSF match and subtract a tem- 
plate from each individual exposure, thus eliminating any 
source that is stationary and of constant brightness. For 
each individual CCD, the template image was the median 
combination of 10 of the best-seeing images. Saturated 
stars were masked at this point to avoid spurious detec- 
tions due to imperfect subtraction. 

To optimally combine the 148 images we defined a 
"weight" for each image based on both its sky back- 
ground and seeing: Wi = cr~^9~^{J2j '^j'^^J^)~^^ where 
9 is the FWHM and a is the standard deviation of the 
background in counts. We defined a reference time to 
be the weighted average of the exposure mid-times. For 
each of the 736 combinations of parallel and perpendicu- 
lar rates we shifted all subtracted images to the reference 
time and computed the weighted average of them. 

We searched each of the "shifted and added" im- 
ages for point sources that corresponded to real or im- 
planted objects moving at the corresponding rate. For 
this we used a combination of the SExtractor pack- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



age (|Bertin fc Arnouti 11996) and a wavelet source de- 
tection routine ( Petit et al.l 200 J ), each with a 3cr de- 
tection threshold. These two approaches rely on very 
different image properties and, thus, have different noise 
characteristics. We have found that the intersection of 
the results from these two routines strongly discriminates 
against false detections. Nevertheless, we find ~ 300 de- 
tections per coadded CCD at each combination of rates 
parallel and perpendicular to the ecliptic. 

We comb through these with an automated search 
algorithm. The algorithm considers all detections in 
all of the parallel and perpendicular rate combinations 
(roughly 2.2 x 10^ detections per CCD), projects them 
on the sky plane, and then searches for clusters of detec- 
tions. Assuming each cluster corresponds to a moving 
object detected in one or more of the coadded images 
for rates close to the real, or synthetic, one, the code 
selects the single detection from every cluster that has 
the smallest ellipticity and largest flux. We visually in- 
spect the resulting ^ 300 candidates per CCD in order 
to reject any obvious spurious detections; those rejected 
are typically poorly subtracted stars that were not ade- 
quately masked, or cosmic rays that resulted in multiple 
detections. This inspection process is extremely fast and 
decreases the number of possible moving objects to ~ 60 
per CCD. We further examine each of these remaining 
candidates to determine the parallel and perpendicular 
rates that yield the PSF that best resembles a normal 
one. In this way we visually found the best rate of mo- 
tion. In this process another 40% of the remaining ob- 
jects were rejected. 

The novelty in this method is the initial automated 
search and visual filter for candidates that are later 
checked by blinking through different rates. The al- 
ternatives are a fully automated search algorithm run 
at a higher detection thres hold to eliminate false posi- 
tives (jBernstein et al.l l2004l ) or a visual search through 
the w hole field and all possible rates (jCladman et al.l 
1200 It) . We expect that a fully visual search would be 
more sensitive, leaving aside the issue of human fatigue. 
To test this we completed a visual search of 20% of 
the area (two of the ten CCDs). In this test we visu- 
ally confirmed or rejected potential candidates identified 
by the two so u rce de tection algorithms, as was done by 
iHolman et al] (|2004D . The comparison showed our ap- 
proach to be ~ 10% less sensitive compared to a visual 
search, roughly independent of magnitude and rate. 

4. CONTROL POPULATION AND DETECTION 
EFFICIENCY 

To calibrate our search we used a control population 
of objects that resemble real TNOs both astromctrically 
and photometrically. We considered a range of distances 
(30 - 200 AU) and magnitudes (i? = 24.5 - 29.5), as 
well as the full range of eccentricities and inclinations to 
create a population of synthetic TNOs. Ephemerides for 
this population wer e created using a modified vers ion of 
the Or6/if routines ijBernstein fc Khushalanill2000( l. 

For each CCD ~ 10 bright stars were used to determine 
a PSF model. These stars are also used to account for 
changes in the seeing and the atmospheric transparency 
when implanting objects. 

We recovered 312 synthetic TNOs and used them to 
determine the efficiency as a function of R magnitude. 
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The result is shown in Figure [T] The detection efficiency 
is well represented by the following function: 



A 
2" 



1 — tanh 



R-Ri 
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(1) 



where the best fit values are A=0.86±0.07, i?5o=26.76± 
0.06 and w=0.38 ± 0.04. The maximum efficiency of our 
search is 86%, and it reaches half this value at magnitude 
R = 26.76. 
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Fig. 1. — Detection efficiency as a function of magnitude. The 
solid line corresponds to Eq. [T] with best fit values: j4=0.86 ± 0.07, 
i?50=26.76 ± 0.06 and w=0.38 ± 0.04. 



The USNO-B Catalog was used to tie the photome- 
try of our objects to a standard system. A number of 
USNO-B stars are unsaturated in our field (m/j > 18). 
The photometric u ncertainty of the catalog is 0.3 mag 
(jMonet et al.ll2003l ). Seven isolated stars were selected, 
and compared to the photometry of an exposure from a 
photometric portion of the night. The flux for a 5- 
pixel (I'.'Ol) aperture over a time t in seconds was found 
to be: 

i? = 27.45 - 2. 5 log /s A. (2) 

An aperture correction of 0.076 mag is included in this 
formula. 

5. RESULTS AND ANALYSIS 

Our search resulted in 20 new objects, summarized in 
Table [TJ The distances to these objects are estimated 
using their parallactic motion, assuming circular orbits. 
To estimate the uncertainty in these distance estimates, 
w e compare the precision of the parallactic rate to that 
in lFuentes fc HolmanI ()2008D and adopt the uncertainty 
in those distances (2.5 AU). There are three objects 
with distances closer than 30 AU; given the poor orbital 
information we assume these objects are Plutinos and 
retain them in our sample. 

We were able to measure the error in our photome- 
try via comparison of the implanted and recovered mag- 
nitudes, plotted in Figure [2l This allowed us to set 
the aperture correction and include any uncertainty in- 
troduced by our method. The net uncertainty is ~0.2 
mag, similar to the expected variance given the sky back- 
ground; we adopted it as the 1-a uncertainty for all our 
detections. As for the accuracy in the rate of motion par- 
allel and perpendicular to the ecliptic is ~ O'.'l h~^. The 
corresponding uncertainty in the angle to the ecliptic is 
~ 3 deg. 

With the standard photometry and efhciency function 
we construct a luminosity function for the 20 objects in 



TABLE 1 

TNO PROPERTIES^ 





Rmag 


R.A. 


Decl. dR.A./dt 


dDecl./dt 


dpar 








(" 


h-i) 


(" h-') 


(AU) 


sdl 


26.7 


22:40:46.967 


-07:28:56.80 


-3.19 


-1.53 


35.5 


sd2 


24.0 


22:40:38.889 


-07:34:33.23 


-2.99 


-1.23 


39.2 


sd3 


24.6 


22:40:48.201 


-07:33:30.04 


-2.99 


-1.23 


39.2 


sd4 


26.7 


22:42:46.043 


-07:25:05.09 


-4.39 


-1.27 


26.8 


sd5 


26.6 


22:42:29.147 


-07:34:07.96 


-2.80 


-1.15 


42.1 


sd6 


25.9 


22:42:39.499 


-07:34:09.84 


-3.55 


-1.46 


32.4 


sd7 


25.4 


22:42:26.506 


-07:34:58.77 


-2.67 


-1.21 


43.6 


sd8 


25.1 


22:42:31.162 


-07:36:21.54 


-2.71 


-1.12 


43.6 


sd9 


25.2 


22:42:36.655 


-07:28:27.49 


-2.28 


-0.83 


53.6 


sdlO 


25.9 


22:42:17.085 


-07:24:26.88 


-3.19 


-1.53 


35.5 


sdll 


24.6 


22:42:18.460 


-07:35:09.91 


-2.89 


-1.19 


40.6 


sdl2 


26.1 


22:41:33.852 


-07:39:08.51 


-2.99 


-1.23 


39.2 


sdl3 


26.6 


22:41:33.660 


-07:40:31.79 


-4.03 


-1.88 


27.5 


sdl4 


26.8 


22:40:57.890 


-07:44:53.89 


-2.89 


-1.19 


40.6 


sdl5 


26.6 


22:41:22.837 


-07:42:07.71 


-4.22 


-1.41 


27.5 


sdl6 


25.7 


22:41:01.823 


-07:35:49.27 


-2.86 


-1.29 


40.6 


sdl7 


26.6 


22:41:09.107 


-07:25:13.78 


-2.78 


-0.93 


43.6 


sdl8 


24.9 


22:41:09.797 


-07:26:31.27 


-2.71 


-1.12 


43.6 


sdl9 


26.7 


22:41:44.583 


-07:31:58.21 


-4.71 


-1.29 


24.9 


sd20 


26.8 


22:42:38.548 


-07:47:45.67 


-2.44 


-1.23 


47.2 



^ Positions and rates valid for MJD 52519.451992. 
found in the text. 



The uncertainties can be 
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Fig. 2. — Histogram of the magnitude difference between im- 
planted and measured magnitudes AR. The dashed line is a Gaus- 
sian of width 0.18 mag. 

this field. The result for this survey is shown on the 
left panel of Figure H The best model for the TNO- 
number cumulative function is given by the harmonic 
mean of two p ower laws, or do u ble po wer law (DPL), as 
int roduced by [B ernstein et all (|2004l ) and corroborated 
bv lFuentes fc H olman (2008). The surface number den- 
sity for the DPL is given by: 

a{R)^C [l0""i(^"^3) + 10(°2-ai)(fl.= ,-23)-a2(fl-23) 

C = S23(l + 10("^-"i)(-""'-23)), (3) 

where ai and a2 are the exponents for the bright and 
faint power law behavior of the model; E23 is the num- 
ber of objects expected brighter than i? = 23 and R^q 
corresponds to the magnitude at which both power law 
behaviors meet. Figure [3] shows the cumulative surface 

density S(i?) = J^^a{x)dx. 

There are not enough objects in this survey alone to 
constrain the model; however the d ata follow the pre- 
vious best estimate for the model (jFuentes fc HolmanI 
l2008f) . The previous best parameters were ai=0.7^Q'^, 
a2=0.3t!];^, cr23=2.0j:gj and i?eg=24.3lgj, overplotted 
in Figure [3| Furthermore, this model predicts 21 detec- 




Fig. 3. — The c umulative number density for the obj ects in our su rvey is shown on the left panel. The solid line is the best previous 
model l IFuentes fc Holman 2008). All surveys listed in Fu entes &: Holm an (2008, Table 2) plus those included in this work are shown in the 
right panel. The most likely model is shown as a solid line (see Fig.|4]l. The function ioO-^3{-R— 23.3) jg ^j^g dot-dashed line. The top axis 
assumes a 4% albedo and every object at 42 AU. The apparent steep slope at i? ~ 26.7 is a result of small number statistics and plotting 
the cumulative, rather thatn differential, luminosity function. 

tions for this search, consistent with the 20 that were 
found. 

We combined o ur s urvey with those hsted in 
iFuentes fc Holmad (|200 8- Table 2). We only considered 
objects found at magnitudes at which the search was 
over 15% of the maximum efficiency. All the objects 
in our search fulfill this requirement. The total lumi- 
nosity function is shown in the right panel of Figure [31 
The DPL best fit is determined throu gh a Markov Chain 
Mont e Carlo simulation (details in iFuentes fc Holmanl 
|2008[) with 10^ steps, and an acceptance rate of 25%. 
The posterior distribution function (PDF) for each pa- 
rameter is plotted in FigureHl the most likely parameters 
are given by: ai=0.73^o°l, a2=0.20lj;;}2, a23=lA6+°-^^ 

and Reg^25.0tos. 
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6. CONCLUSIONS 

We have performed a pencil-beam search of a single 
Subaru ecliptic field. The search covered 0.255 deg^, with 
a limiting magnitude R = 26.76. We found 20 new TNOs 
wi th magnitudes between R = 24.0 and 26.8. As argued 
by iBernstein et al.l (|2004f ) , it is not surprising that our 
faintest detection is near our 50% efficiency threshold, 
given that the luminosity function is much shallower at 
R - 27. 

Including other surveys in the analysis (see § [5]) we de- 
rive the most likely DPL model for the luminosity func- 
tion. With only 20 new objects there is an important im- 
provem ent in the model's PDF from IFuentes fc Holmanl 
()2008l Figure 13) to the current result (See Figure g]); 
due to the new detections constraining the model be- 
tween i? = 26 and 27. 

The bright end power-law exponent Q;i=0.73tQ 09 is 
very close to the re s ults of previous shallower sur- 
veys (|Gladman et all 1200 If ). The break magnitude 
i?eg=25.0lQg is more than 1-a fainter than the ini- 
tial estimate by IBernstein et alj ()2004l ). This is also 
consistent with earlier surveys reporting excellent fits 




Fig. 4. — The PDF for each model's parameter. The most 
likely value for each parameter and 68% confidence limits are: 
ai=0.73+°-°l, a2=0.20to\l, '^23=i-i6l:o\2 and fie9=25.0^ 
These are plotted as vertical lines. 
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to a single power model for magnitudes brighter than 
i? ^ 26 (Gladman et al. 2001). Further data between 
the current ground-based detection limit of i? 27 and 
the HST's at i? ~ 28.5 will determine the break magni- 
tude even more accurately. 

Assuming that all objects are located at an heliocentric 
distance of 42 AU, the break in the luminosity function 
Req reflects a break in the size distribution D. The corre- 
sponding diameter at which the size distribution breaks 
is £» = 90 ± 30(p/0.04)-°-5 km, where p is the albedo. 
Theo ries predict the existence of a break (iKenvon et al.l 
12008( 1 ■ but at smaller diameters (D ~ 20 - 40 km). The 
difference reflects the uncertainty in the initial conditions 
for numerical simulations and the observational assump- 
tions regarding physical properties like the albedo. 
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In this survey both classical and excited TNO popu- A comparison with other populations in the solar system 
lations are entangled. With better orbital information can shed light on their origin, 
the size distribution of each can be studied separately. 
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